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CIRCUMFERENTIAL INLET DISTORTIONS IN AN AXIAL FLOW COMPRESSOR 



by 

Fraiik R. Carter 

Submitted to the Department of Aeronautical Engineering on May 26 , 
1958^ ill partial fiLLfillment of the requirements for the degree of Master of 
Science. 



ABSTRACT 



Due to the importance of inlet velocity profile distortions in 
axial compressors of high-speed, gas txirbine driven aircraft, this study 
haa been made to investigate velocity inlet defects and their attentuation 
for the case of an isolated rotor in an axial flow ccm^iressor. The effects 
of various profile generators at several flow rates were investigated, and 
the actual measvired effect of the rotor on the inlet velocity distortion 
was compared with the theoretically predicted downstream profile. The theory 
used was an incompressible, invi sc id actuator disc analysis with the assvunp- 
tion that the hub/tip ratio is close to \mity, reducing the problem to one 
of two dimensions. It was found that the velocity defects were attentiated 
across the rotor since the distorted flow is energized more than the undis- 
torted flow by the rotor. Further, it was found that the theory proposed 
by Ehrich did predict the mean streamline shift with reasonable accxiracy as 
well as the general shape of the downstream profile within the limits of 
the calculations and tests performed. Refinements in the theory, however, 
do appear to be in order, as the theory undexestimated the distortion leav- 
ing the rotor by as much as 25 fo. 
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SYMBOLS 



Symbols 



a 


tangent a 


b 


tangent p 




eucial gas velocity 


D 


mean diameter of compressor 


P 


static pressiire 




total pressure 


Q 


volume flow rate, ft^/ sec 


q. 


dynamic pressure 


u 


wheel speed 


a 


absolute flow angle 


P 


relative flow angle 


P 


Cjj/U, mean flow coefficient 


0 


circumferential phase orientation of distortion 


Q 


angle measurement in circvanferential direction 


P 


yaw angle 



Subscripts 



1 


some distance upstream of the blade row 


2 


immediately in front of the blade row 


3 


immediately downstream of the blade row 


h 


some distance downstream of the blade 



1. Introduction 



Inlet velocity profile distortions in axial compressors have 
been a most important problem in high-speed gas tvirbine driven aircraft 
as these distortions reduce the performance of the over-all power plant 
by reducing the performance of the compressor. Such distortions may oc- 
cxu: in flight under various conditions such as at off- design conditions, 
when the aircraft is flown at high angles of attack, or when the flow is 
turned or diffused too rapidly within the diffuser causing wall separation 
within the diffuser ducting. A distorted flow may also res\ilt because of 
a non-uniform compression at the inlet throat at supersonic, speeds when 
the engine is operated at off- design conditions or when the aircraft is 
flown at high angles of attack. Further, circxmiferential inflow distor- 
tions have proved to be one of the most significant causes for the re- 
duction in the stall margin of an axial compressor, and thus the study of 
such inlet conditions, their prevention, reduction, and attenuation has 
been the subject of several recent analytical publications, though little 
work has been done experimental 1 y . 

In the design of a compressor an understanding of the effects 
of a distortion and its attentuation makes it possible to design a compres- 
sor which is better able to withstand these effects. An investigation of 
the velocity diagrams of representative inlet stages shows that, in gen- 
eral, velocity profile defects will be reduced across tlie rotor because 
the distorted flow, or low energy region, is energized more than the un- 
distorted flow by the rotor. The pxirpose of this present experimental 
investigation is to study the case of a single isolated rotor in an axial 
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flow compressor at different flow rates and with different inlet velocity 
profiles. The actual effect of the rotor on the inlet velocity profile 
was measiired and the result compared with the do\mstream profile predicted 
hy theories set forth in papers by Fredric Ehrich of the General Electric 
Company and George C. Ashby, Jr., of the NACA. 

Erich's approach, which appears to be particularly satisfactory, 
is an antuator disc analysis. The self transport of vorticity is neglected 
in order to linearize the equations, and the analysis is restricted to the 
study of incompressible, inviscid flow with the assumption of a hub/tip 
ratio close to unity reducing the problem to two- dimensions. He had de- 
rived expressions giving the attenuation of the velocity distortion through 
the rotor and the circumferential shift in the profile as a function of the 
rotor geometry. 
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2 . Equipment 



2.1 Compressor 

The compressor used in these tests and the test set up are shown 
in Figure 1. A single isolated rotor was used in all tests, and the char- 
acteristics of the compressor are presented in Reference 1. The compressor 
geometry is; 



Outside radius 
Huh/tip ratio 
Blade chord 
Blade aspect ratio 



U ..625 inches 

0.75 

1.00 inch 
2.84 



The blade-tip clearance was approximately 0.055 inches, and the blade sec- 
tions were NACA 65- (12) 10 £U.rfoils with a stagger at the mean radius of 
52.7 degrees measured from the axial direction. The cascade solidity was 

1.02 at the leading edge, and there was linear twist of 9*70 degrees from 
the root to the tip. A schematic diagram of the compressor with the avail- 
able press\ire reading stations is shown as Figure 2. 

2.2 Profile Generators 

I'o generate the various inlet velocity profiles, screens of var- 
ious solidities were used. Tests were first conducted* with 22 1/2, 45, 
and 90 degree screens of high solidity, shoim at the left of Figure where- 
as the majority of the tests were conducted with one or two layers of window 
screening in a 45 degree screen holder, as shown on the right of Figure 5 
with the screen holder, ctirved fences, and extension arm. 
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2 . 3 Instrument at ion 

The probes used for measiaring yaw angle, static and total pres- 
sures were of the three hole Fecheimer tube type described in Reference 2. 
These probes haVe three 0.012 inch holes drilled through th^ roimded end 
of the tube at spacings of 42 degrees of arc with three hypodermic tubes 
contained inside a 0.125 inch stock tube. The middle hole points into the 
flow when the pressures at the outer two drillings are equalized. Probe 
calibration procedures and results are contained in Appendix A. 

The yawing of the probes was accomplished through the use of a 
standard strain gage pressure transducer manxtfactured by the Statham Lab- 
oratories, Inc., of Los Angeles, California. Stathpm pressure transducers 
offer a simple means for the conversion of dynamic or static pressure into 
sensible electrical units. They are accurate and linear within ^ 1^ of 
fvill scale. The ya|w angles were read with protractors mounted as shown in 
Figure 4 on the compressor shell. These protractors were aligned with a 
square uping a known axial, flat surface on the outer wall of the compres- 
sor as a reference. The probes were aligned using the following optical 
technique. A small mirror was mounted on the, shaft of each probe as seen 
in Figure 4 with the plane of the mirror being perpendicular to the axijp 
of the center hole. A reference mirror was then aligned with the center 
line bf the compressor using the axial, flat surface mentioned above. The 
mirror had a scribed vertical line and a "peep" hole which were used to 
bring into coincidence all images of the scribed line when looking throiagh 
the "peep" hole at the mirror mounted on the probe. With this accomplished 
it was assured that the probe was aligned 90 degrees to the centerline of 
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the compressor and coiold then he locked into position with respect to the 
protractor. Similar procedures are employed in aligning interferometer 
mirrors. 

Pressxure measurements upstream and downstream were mads throu^ 
the use of a manifold pressure box, shown in Figure 1, and a strain gage 
pressxure transducer connected to a D-C calibrating device which incorpo- 
rates a transistorized D-C amplifier. This transducer equipment was manu- 
factxured by the Dynamic Instirument Company of Cambridge, Massachusetts. 
Some of the tests were rxm using a standard strain gage pressure trans- 
ducer manxifactxured by the Statham Laboratories. A 32- tube inclined ma- 
nometer board having a slope of 7*6 degrees and using red oil (specific 
gravity, O. 827 ) 'was used as stand-by equipment and for checking the pres- 
sure equipment described when necessary. This board was calibrated with 
a Prandtl micromanometer. 

GSie damper position and the r.p.m. of the rotor controlled the 
volxone flow of air throxxgh the con^jressor, and they were measured through 
the use of a coxinter moimted on the compressor control panel and an elec- 
trical strobotac respectively. 



It' 




3. Investigation of an Inlet Total Pressure Distortion of Sma3.1. 
Circumferential Extent 



3«1 Experimental Procedure 

For the experimental tests, the total pressvire distortion gen- 
erated by a one- inch diameter rod was measured upstream and downstream of 
the rotor. By using the measured upstream profile and a derived analyti- 
cal equation for incompressible flow derived by George C. Ashby, Jr., in 
Reference h, the profile downstream of the rotor was calculated and com- 
pared with the measured downstream conditions. The one- inch diameter rod 
was installed four inches upstream of the blade row of the rotor, the vip- 
stream probe one inch upstream of the rotor, and the downstream probe one 
inch downstream of the rotor. The cylinder was moved circumferentially 
throu^ the use of an electrical traversing rig while the probes ronained 
fixed in position. Readings were taJcen at two degree intervals of traverse. 

3.2 Results and Discussion 

Equation (2) of Reference 4 derived for incompressible flow is 

cos ((X, + ) cos + Bp) cos aicos Bg 

cos Bicos Ctp 

Ashby states that when the total pressure deficit between the undisturbed 
flow and a point in the disturbed flow upstream of the rotor and the ve- 
locity diagram for the xmdisturbed flow are known, the doimstream total 
pressure deficit can be determined using the above equation. The distor- 
tion must be of small circxanferential extent however, so that the im- 
distorted flow dictates the static pressure field in the blade passage. 



^Po2 = ^ol 
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Figure 5 shows the measixred upstream total pressure profile and 
the measixred and predicted profiles downstream of the rotor. This theory 
failed to predict the downstream profile with any reasonable accuracy. Fig- 
ure 4 of Reference 4 showed similar resvilts with the measured downstretam 
profile showing a larger attenuation and wider circumferential disturbance 
than was predicted by the theory. 

Ashby asserts that reasonable accuracy is obtained when the static 
pressure field in the blade passage is altered by a distortion of large cir- 
cumferential extent. That this is not the case has been determined by sev- 
eral tests, one of which is exemplified by Figure 6 which shows the meas- 
ured total pressixre distortions upstream and downstream of the rotor and 
the calculated downstream profile using a 45 degree window screen. 

$.5 Conclusions 

A single experimental investigation of the effect of an isolated 
rotor on a known inlet total pressure distortion in an axiel flow compres- 
sor and a cotDparison of this effect with that predicted by a theoretical 
equation derived in Reference 4 has been conducted. The predicted down- 
stream conditions compared unfavorably with the measured conditions, and 
therefore it was concluded that Ashby's theory as set forth in Reference 
4 is not practical for predicting significant inlet distortions. 
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k. Velocity Profile Generation 

In the first experiments conducted, the screens of various sol- 
idities were tested. These screens were traversed manually with the probes 
left stationary. Upon taking data however with the 45 degree screen of 
highest solidity, it was found that the direction of the flow upstream of 
the rotor was changed radically from essentially axial to a maximum of 
sixty degrees from axial. The restilting inlet velocity profile is shown 
on Figure 7 . It was noted that the defect was not symmetrical, having a 
larger velocity on the side down from the direction of rotation. The 
cause and remedy of this asymmetry was later invest igaljed. It was also 
noted in this profile that there were two larger reductions in the axial 
velocity at the edges of the distortion giving the profile an appearance 
of a "W". As the screen holder width is one-eighth of an inch and is 
located three to five inches upstream frcan the measuring probe, it was 
concluded that the W-sh§pe was not caused by the sides of thS' screen 
holder as any wake caused by the holder would be dissipated prior to 
progressing downstream to the probe, 24 or more chord lengths away. As 
the tests with and without fences both give the characteristic W- shape, 
it was further concluded that the curved fences which may cause thicken- 
ing of the bovundary layer did not cause this iibreguiLar shape. As the 
width of the distortion was roughly the width of the screen, it was de- 
duced that the edges of the screen (or the screen holder) were generat- 
ing vortex sheets which were extending downstream causing the increased 
defect due to the direction of the circulation in the vortex sheet. This 
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is in agreement with theory as stated in Reference 7 which shows a plot 
of velocity versus distance for a pair of point vortices as 




Tests were discontinued with these screens, as the profiles generated would 
be unsxiitable for comparison with a l pertxarbation theory used in Ref- 
erence 5 to the large divergence of the flow. 

To reduce the blockage and consequently the absolute angle of 
flow, window screening was used. Again, however, the aiigles were too large 
to be considered within the limits of small perturbation theory. These lim- 
its were set as five degrees either side of axial. 

Fences were attached to the screen holder in a manner shown in 
Figure 5 to reduce the divergence angle. This broxi^t the angles within 
ten degrees of axial, and finally less than five degrees when the fences 
were c\irved upstream to align them roughly with the streamlines to pre- 
vent the divergence of the flow at the leading edges of the fences. The 
fence was first attached with twelve inches extending upstream and noth- 
ing downstream giving the profile shown on Figure 8. The straight fence 
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with twelve inches extending upstream and six inches downstream gave the 
resulting profile plotted on Figure 9» 

A further modification which in5)roved the results was the shorten- 
ing of the upstream fence to five inches and extending the downstream one 
inch. With this design;, several tests were made for one ^d two layers 
of window screening, various flow coefficients, and volume rates. These 
results are shown on Figures 10, 11, 12, and I5. These figures all showed 
the asymmetry of the first profiles and the characteristic W- shape due to 
the vortex sheets emanating from the edges of the screen. The probe loca- 
tions, the volvone flow rates, the flow coefficients, and certain remarks 
for all the runs are contained in Table I for convenient reference with 
Figure 15 explaii^ing colimm headings. 

Due to the proximity of the rotor to the Screen and probe, it 
can be concluded that the asymmetry of the various profiles was due to 
rotor effect. Since actuator disc theory shows that rotor influence di- 
minishes exponentially with distance from the rotor, an att^pt was made 
to reduce the asymmetry by extending the screen upstream eight inches and 
moving the probe upstream ten inches. The screen was extended by attach- 
ing an arm shown in Figure 5 of length eight inches. The asymmetry was 
slightly reduced, but not by a significant amount. A further attempt to 
solve this problem was the installation of two screens I80 degrees apart 
from one another at this extended position. The theoretical effect of 
this attempt is contained in Appendix B. The profile generated by this 
arrangement is shown on Figure l4, and shows again that thi rotor effect 
is still present. It was then concluded that a symmetrical profile must 
be generated by modifying the screening to compensate for the effect of 
the rotor rather than trying to eliminate it. 
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3. Investigation of Velocity Profile Distortions of Significant 
Circvanferential Extent 



$♦1 Analytical Procedure 

From the several inlet distortions of Figures 7 through 15 ^ two 
were selected to be analyzed using the theory proposed by Ehrich in Refer- 
ence 5* These two were those shown in Figures 11 and. I 5 with flow co- 
efficients of 0.440 and 0.558 respectively. In the test of Figure 11 two 
layers of window screening were used and in that of Figure I 5 , one layer. 

With the W- shape inlet profiles obtained by tests, it was found 
that more Fourier coefficients were needed to duplicate the inlet profile 
thian was deemed practicable to handle in the analysis. Therefore, as this 
W- shape was caused by vortex sheets which probably wash out or which are 
undoubtedly washed out by the rotor, an assvimed smooth curve was faired 
in thus eliminating these irregularities caused by the vortex sheets. The 
assmed alterations of the inlet profiles are shown as dashed lines on Fig- 
ures 16 and 17. 

The Fourier coefficients were obtained by using the Dent-Draper 
rolling sphere harmonic analyzer manufactured by the Mico Instrument Co. 
of Cambridge, Massachusetts. It was found that the inlet profile could 
be duplicated closely at most points with ten harmonics, but that sev- 
eral of the points required twenty harmonics bo reproduce the c\arve with 
reasonable accxiracy. These calculated points are shown on Figures I 6 and 
17 as crosses, and it can be seen that they reproduce the profile with 
about ten per cent accuracy. It was assumed that this accuracy was main- 
tained in the analysis for finding the downstream profiles. This seemed 
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reasonable as the inlet coefficients were only operated on by analytical 
expressions set forth in Reference 5« 

Figures l6 and I7 show the measured upstream and downstream pro- 
files, the calculated downstream profile, the experimental points for the 
upstream and downstream profiles, and the calculated upstream and down- 
stream points used in this analysis. 

$.2 Disdussion 

In Ehrich's theory, an ejq>ression for tan 04 is given. This ex- 
pression, however, does not allow one to determine the principal values of 
04« .Small, changes in the denominator of the expression for tan 04 when 04 
is near ^90 degrees can cause the sign of the tangent to cha4ge. As small 
changes should not cause large phase angle differences such as from +90 de- 
grees to -89 degrees, it was concluded that the angle must either lie id.thin 
the first and second quadrants or within the third and foxurth quadrants, 
either region allowing tan 04 to take on all values from - 00 to + co . As 
Reference 5 gives an example of -40 degrees for the phase angle, it was 
concluded that Ehrich meant for the angle 04 to range only from 0 to -I80 
degrees. This is also in agreement with one's intuition in regards to phase 
shift in a three-dimensional case with a rotor of finite thickness. 

From Figures I6 and I7, it is apparent that the theory does pre- 
dict the phase shift with a considerable degree of accuracy. The general 
shape of the profile also appears to be correct but the magnitudes were not 
predicted with similar accxnracy. Part of this difference may lie in the 



fact that too few harmonics were used in the calculations. 
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3»5 Conclusions and Sijggestions for Futvire Work 

It is believed that Ehrich's theory predicts with reasonable ac- 
curacy the phase shift of the mean streamlines throu^ an isolated rotor. 

It further appears that his theory predicts the general shape of the down- 
stream profile, but may need certain refinements to bring the magnitude 
of the defect into closer agreement with the measured defect. This may 
be rectified in futvire work using higher flow coefficients and varying 
solidity screens, as it was noted that in Figure 17 with a flow coefficient 
of 0.558> closer correlation resulted than in Figure I6 with a flow coef- 
ficient of 0.4k0. In the example used by Ehrich, the flow coefficient was 
equal to 1.00 indicating that the tests conducted may be near the limits 
of ^ for the theoretical predictions to be valid. 

It is therefore suggested that future tests be conducted with 
the results of this thesis as a basis for the selection of flow coeffi- 
cients, volume flow rates, rotor speeds, screening to be used, and loca- 
tions for taking measurements. It is further suggested that the Fourier 
coefficients be calculated on the IBM 70^ computer (for example, E2PK 
4EA1, Fourier series coefficient sub-routine) to increase the accuracy 
of the computations and to reduce the laborious calculations required in 
operating on ten to twenty harmonics. To help reduce the number of terms 
required in the seried, it is suggested that the inlet distortion be al- 
tered by making a screen of varying solidity which would eliminate the W- 
shape and perhaps make the profile symmetrical. 
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table I 



Fig. 


Obstruction 


a 


b 


C 


d 




Q sec 


Remarks 




















10 


45* window screen 
(two layers) 


3" 


14" 


9" 


390 

arc 


0.319 


36.7 


5 "curved fence up- 
stream and 1" cxurved 
fence downstream 


8 


45°window screen 
(one layer) 


5" 


4" 


9" 


0“ 


0.364 


36.6 


12"curved fence up- 
stream and no fence 
downstream 


9 


45 “window screen 
(one layer) 


5" 


4" 


9" 


0“ 


0.391 


39.4 


12"straight fence up- 
stream and 6" fence 
downstream 


15 


Two 45 “window 
screens separat- 
ed by 180“ (one 
layer each) 


3" 


14" 


9" 


390 


o.4l^ 


48.8 


No fences 


11 


45 “window screen 
(two layers) 


3" 


14" 


9" 


39" 


0.440 


50.6 


Same as Fig. 5 


12 


45 “window screen 
.. (two layers) 


3" 


14" 


9" 


39" 


0.48l 


48.4 


Same as Fig. 5 


7- 


45“ screen 
(highest solidity) 


5" 


4- 

1 


9" 


0“ 


0.501 


60.5 


No fences 


13 


45 “window screen 
(one layer) 


3" 


14" 


9" 


39" 


0.538 


54.3 


Same as Fig. 5 
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APPENDIX A 



PROBE CALIBRATION 



Probe calibration was carried out in the Gas Turbine Labora- 



tory probe calibration tunnel which was used in the experimental work of 
Reference 6. Figure l8 shows the total pressure gradient measvtred. This 
measured profile was in close agreement with the profile obtained from Ref- 
erence 6, as can be seen on Figure l8. 



until the static pressvire taps were equalized. The resiilting yaw angle, Kjj , 
the total pressxire, and the static pressure were read, with the pressures 
beiiag read on the Dynisco pressure reading apparatus. 



obtained from Reference 6, was determined to be and 5 *72 for probe 

"B". The angles obtained using the three-hole probes checked extremely 



closely with angles obtained in the experiments of Reference 6. As the 
total pressure gradient was measured in the compressor and was found to 
be small, (relative to the size of the probe) a correction for the yaw 
angle measurements was tmnecessary. Also, since the total pressvire gra- 
dient in the compressor was small and the yaw angles ancvirate, corrections 



The procedure used in measurement was to yaw the three-hole probe 



The constant for probe "A" in the e:q>ression for yaw angle cor- 




/ 
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for total pressiure were considered negligible. The static pressvires, 
however, were corrected, as the static holes in the Fecheimer tube were 
not drilled at the location at which free stream static pressxire theoreti- 
cally exists. Machining problems prevented the choice of this location. 
Figure 19 shows a plot of (Pindicated'^wall/*^-^ versus vertical traverse 
in the calibration timnel. This plot gave the results that for probe "A" 

^true ~ Pprobe = O.I52 

For probe "B": 

Ptrue "Pprobe ^ 0.O85 
<1 






! 
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APPENDIX B 



CHOICE OF DISTORTION WAVE LENGTH AND ITS DISTANCE FROM THE ROTOR 



In the analytical investigation of Reference 5 > it was shown 

that any circxanferential distortion of width, 2 jtr/n, will induce local 

distxarhances at the rotor. These disturbances will die away to a value 

of exp - 2 jt upstream and downstream of the rotor in an axial distance 

equal to the width of the distortion . 

For the tests conducted the distortion width is 

distortion width = 2 jtr/n = 20 jt 

therefore r _ 
n 

The decay term in the attenuating process is 

_ nx 
e'^ 

and thus with the probe located four inches upstream of the rotor, the 
rotor influenpe, or disturbance, is reduced by 

exp (-O.h) = 0.67 

Increasing x to nine inches by moving the probe upstream causes the dis- 
turbance to diminish to 

exp (-0.9) = O.U06 

An alternative method for decreasing the rotor influence is to 
change the period or distortion width by introducing a second screen Id- 
eated 180 degrees away from the first. Thus: = lOjt £ = 5 

and with x = 9" the distvirbance diminishes theoretically to 

exp (-1.8) = 0.1651 
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APPEHDIX C 

A SUMMARY OF EHRICH'S THEORY FOR AN ISOLATED ROTOR 



From Reference 5 ^ Tor an inlet profile defined by 

Cx, - G ^ 2(^.U -(nj 

the velocity at the exit for an isolated rotor (referred to the same 
mean streamline) will be given by 

where the rotor attenuation is given by 

\(:J J 

and the phase angle shift is given by 
vdiere 




ajL = tan 

bg = u/c^ 

(bj - aij) = U/o'x 

For the data of Figure l6; 

9= 



£14 = tan 

bj = bg + ai^ - a^ 



For the data of Figure I7; 




0.7077 



^4 = -75"6 



= - 88?9 
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FIG. 1 — SET-UP OF TEST APPARATUS 




FIG. 2 — DIAGRAM OF THE COMPRESSOR 
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FIG. 3 — PROFILE GENERATORS 
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FIG. 4 — PROBE ARRANGEMENT 
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FIG. 5 “ TOTAL PRESSURE PROFILE (CYLINDER) 
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FIG. 8 — UPSTREAM PROFILE 
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FIG. 10 — UPSTREAM PROFILE 
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FIG. 12 — UPSTREAM PROFILE 





FIG, 14 — UPSTREAM PROFILE 
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FIG. 15 — MEASUREMENT REFERENCE POSITIONS 
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FIG. 16 — MEASURED AND PREDICTED PROFILES 
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FIG. 17 — MEASURED AND PREDICTED PROFILES 
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FIG. 18 — CALIBRATION TUNNEL TOTAL PRESSURE GRADIENT 
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FIG. 19 — PROBE CALIBRATION 
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